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ABSTRACT 

The problem of small perturbation potential supersonic flow 
around complex configurations is considered. This problem requires 
the solution of an integral equation relating the values of the 
potential on the surface of the body to the values of the normal 
derivative, which is known from the small perturbation boundary 
conditions. The surface of the body is divided into small (hyper- 
boloidal quadrilateral) surface elements, , which are described 
in terms of the Cartesian components of the four corner points. 

The values of the potential (and its normal derivative) within 
each element is assumed to be constant and equal to its value 
at the centroid of the element. This yields a set of linear 
algebraic equations. The coefficients of the equation are given 
by source and doublet integrals over the surface elements, iL ^ . 

The results obtained using the above formulation are compared 
with existing analytical and experimental results. 
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INTRODUCTI^N 

The evaluation of the aerodynamic loads is important in 
structural design and flight control of the aircraft. Those 
oadn must he evaluated by nn iterative process which requires 
that the mathematical modeling of tne problem bo general, 
flexible and efficient. The usual methods for the evaluation of 
aerodynamic loads is the computational lifting-surface theory, 
these methods are efficient and flexible but not general 
enough to consider problem for complex configurations. Some 
computational methods around complex configurations have already 
developed. Ho-.e/cr, they are usually quite cumbersome to use. 
Furthermore, for oscillatory flow around complex configurations, 
only techniques based on the doublet -? attice method exist in 
subsonic range while no method is available in the supersonic 
ones . 

The development of the present method is aimed to overcome 
the shortages of the existing methods and provide an efficient, 
general and flexible aerodynamic tool to be used in structural 
design and flight control of the aircraft. 

The present method is mainly based on the theoretical 
formulations developed by Morino (Refs. 1, 2, and 3). The 
geometry of the aircraft and wake is approximated by a number 
of quadrilateral elements described by hyperboloidal surfaces. 

By applyinq the Green function method, one obtains a linear 
equation relating the velocity potential , at any point, , 
in the flow field with the values of <j' and its normal deriva- 
tives on the surface /. , surrounding the body and the wake. 
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INTRODUCTION 

Tho evaluation of the torodynamic loads is important in 
structural design and flight control of the aircraft. Those* 
loads must be evaluated by an iterative process which requires 
that the mathematical modeling of tho problem bo general, 
flexible and efficient. The usual methods for the evaluation of 
aerodynamic loads is tho computational lifting-surface theory, 
these methods are efficient and flexible but not general 
enough to consider problem for complex configurations. Some 
computational methods around complex configurations have already 
developed. However, they are usually quite cumbersome to use. 
Furthermore, for oscillatory flow around complex configurations, 
only techniques based on the doublet - 1 attice method exist in 
subsonic range while no method is available in the supersonic 
ones . 

The development of tho present method is aimed to overcome 
the shortages of the existing method}! and provide an efficient, 
general and flexible aerodynamic tool to be used in structural 
design and flight control of tho aircraft. 

The present method is mainly based on the theoretical 
formulations developed by Morino (Refs. 1, 2, and 3). The 
geometry of the aircraft and wake is approximated by a number 
of quadrilateral elements described by hypcrboloidal surfaces. 

By applying the Green function method, one obtains a linear 
equation relating the velocity potential t|’ , at any point, p^ , 
in the flow field with the values of c|' and its normal deriva- 
tives on tho surface /. , surrounding the body and the wake. 



An integral equation is obtained by imposing that the value of 
the potential at approaches the value of that at point p on 
the surface if approach p . Tho value of the veloclt , poten- 
tial, u is assumed to be constant within each element. Then, 
the integral equation becomes a system of algebraic equations 
which relates unknowns, if ^ , at the centroids of elements, 
with coefficients evaluated analytically. Once the distribution 
of the velocity potential is obtained, the pressure distribution 
as well as the generalized forces are evaluated. 

In fection 1, tho basic equations are introduced. In 
Section 2, the numerical formulation is presented. In Section 
3 Jnite values of tho integral are considered. Oscillatory 
..lows are discussed in Section 4. In Section 5, numerical re- 
sults are presented. In Appendices A, B, and C, useful equations 
are derived. 


SECTION I 


BASIC EQUATIONS 

1 . 1 I ntroductio n 

In this section, basic flow equations arc first introduced 
in Subsection 1.2, while boundary condition of the problem is 
considered in Subsection 1.3. The Green theorems for steady 
and oscillatory flow arc considered in Subsection 1.4. The in- 
tegration scheme and the hyperboloidal .. lement are considered 
in Subsection 1.5. In Subsocion 1.6 the role of the diaphrapms 
is discussed. Finally in Subsection 1.7 the evaluation of the 
pressure coefficient and the generalized forces are considered. 


1,2 Basic Flow Equati on 

The flow is assumed to bo isentropic, inviscid and initially 
irrotational such that the flow can be described b} the velocity 
potential ‘I’ . Consider a frame of reference such that the un- 
disturbed flow has a velocity in the direction the the posi- 
tive x axis, so that 

5 x. + <. ( 1 . 1 ) 


whore ij is the perturbation velocity potential. Then, the 
linearized equation of potential flow based on the assumption 
of sir.* 11 ! perturbation is given by 

V Z d> - 1 flli = O (1.2) 

' c\ r cUi L 

. jre 


i , 2 . 

dX 

is the linearized total time derivative. 


. o _L 

+ “ 3* 


(1.3) 


(1,4) 


Consider the nondimoni,xonal variables 

X,</hL Y- tit 

whero 

B /IT**- I ’ U.5) 

and t is a reference length, 

1,3 Boundary Conditi on 

The lifting body considered hero lias complex configura- 
tion and is moving with small vibration with arbitrary mode. 

Thus the surface of the body is represented in the general form 

£ (x,y,?, i) = 0 »•«) 

and the boundary condition on the surface of the body is given 
by 

? l l 4 V<3>. V& z 0 (1.7) 

at 

By using Eq. (1.1), Eq. (1,7) becomes 



Furthermore, because th** flow is uniform at infinity, the boundary 
condition at infinity is given by j) = 0. For steady subsonic 
flow, Eq. (1.8) is simplified as 

. n z - u (1.9) 

where Y\ is the surface normal of the aircraft and n is the 
x - component of , 

In order to use compact vector notation, it is conven- 
ient to introduce the concept of the conormal and cogradient 

* y J 1 


( 1 . 10 ) 


( 1 . 11 ) 


y 4 . K . 1 T - i. k 

“ ax ar J az 

Furthermore, the conormal derivative is given by 

£ . ft c . V =• V e 

i >ij c ~ 

With this notatior. , Eq (1.9) becomes. 


(1.12) 


or 


fk 


(is U + 2S ?^U' 

' a.c ;>x 1 3;' 3/ + 'ii ) * 

V 0 

. u,/. i . 1 . ?£ , I as'i 

' tTv®*»X •>; ;r ay 3* 5z E 35? / 

« Mj ff. £§ ■ ? ;• , ; 1 4 1 \ £§ , f 1 21 

*(. L v >X < i ' ft / • ft J 

. w 1 lvsl = o 

e L*ii‘ & 'Mb j ' ' 


?-t - - n,(i + 23i .. 13) 

3N C ' b it J (1<13 * 

The second term on the right hand side is negligible . Tk us 


r. - (1.14) 

In addition to the above notation, it is convenient to introduce 
a ppecial algeora, called supersonic vector algebra or, super - 
algebra , to simplify the algebraic manipulation for the supersonic 
flow theory. Details of this algebra and the proofs of some 
rules of this algebra are given in Appendix A. 


1,4 Green's Theorems 

According to Eg. (1.2), the Green function is given by 



- 4 - 


-J — * r .,^ Lt ) 


- 77 , ’ w f, »' ,, - , (115) 

*1 Cv\ 

•* 

with 

3 -- O . (1.16) 

at infinity. 

Detail of tha derivation of function G is presented in Ref. 

1. The results are summarized here. Tor steady supersonic flow 
G is given by 
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2r, r 


where H is given by Tq (1,22) while 


„ Mi) '/i 


For oscillatory supersonic flow G is given by 
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4n 
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Wi'ere 


Q ~ - -L r Ai (x- * j + r 
-V.5* L J 


(1.17) 


(1.18) 


(1.19) 


( 1 . 20 ) 


By applying the Green function method and using Eq. (1.17), 
the linearized equation for the steady supersonic flow can be 
derived as 


ivll 


- £ M 4 « * & * ) * 


i 


f; el n U/ii i A C 7l 

r O •>’ f.'dc £ 


H.J for *.-X > 

. 0 f.r x„-X < [(V.-Y)\fo-Z)‘] 

1/ /?// is the supernorm of vector R ’ namely 

ft! 


( 1 . 21 ) 

(1.2 2 a) 

(1 22 t) 
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with "super-produce", O , defined by (see Appendix A) 

A® b ex* L, - A y fcy - 

and the conormal derivative is given by Eq (1,12). 
For op dilatory flow, the Green theorem is given by 

a a j }(K - - <£> il JL 

1 ^ am if R 3 

r & ^ 

A 

where is given by 


, , / ly , _\ -.'Jiir- M/j 

^ - <£(*^,* 2 , TJ e 


(1.24) 


(1.25) 


(1.26) 


1,5 Numerical Procedure 

Dy imposing that the value of the potential at ^ approaches 
the value of potential at p on the surface X ,• if ; , approach 
p t t. e va! ae of £ ( | ) »a Eqs (1.21) and (1.25) is found to be 
1/2 (see App. F ; V.i f 2). Then, an integral equation relating the 
potential on the surface 7 to its normal derivative is ob- 
tained. In order to solve this integral equation the surface, 

/L , surrounding the aircraft and wake and diaphragm, if necessary 
is divided into a number of quadrilateral elements which arc 
approximated by hyperbo] oidal elements. The general expression 
of hyperboloidal elements can be written ai 

p-. p, . \ p .rfp. + s,^p 3 

whore. j\, j>, , ^ , and J' are .linear combinations of the four 


(1.27) 


corner vectors, j^ 4 t „ t , » and p_ . (see fig. 1) 

By assuming that the values of velocity potential and its 
normal derivatives within each element, 1 , are constant, 

Eq. (21) for steady supersonic flow reduces to 



where ( ) is any point on the surface _ , and N is the number 

‘ 

of elements on _ . 

Imposing the condition that Eq. (1.28) is satisfied at the 

centroids, p of the elements, l , a system of algebraic 

i h 

equations ii obtained. This yields 
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(1.29) 


(1.30) 


(1.31) 


and the definition of >J is similar to the one for subsonic 


fc W 

flow (see Ref. 1). However, in this report, the wake is not 
included since only supersonic-trailing-edge configurations 
are considered here. 

Equations for oscillatory supersonic flow are considered 
in Bection 4. 


1 . G Diaphragm s 

There are three categories of win.j geometry in supersonic 


flow. The first one. are wings with supersonic leading edge. For 
Chis kind of wings, the clement, , on the upper (or lower) 
surface are influenced only by the elements on the surface of 
the same side, therefore, the integrations in Equations ( 1 . 30 ) 
and ( 1 . 31 ) are performed over the area o i one side only. The 
second one are wings with subsonic leading edge. For this kind 
of wings, diaphragm may or may not be used. For wing3 of large 
thickness, results obtained with or without diaphragm are al- 
most the same. However, for wings of snail thickness, diaphragms 
arc suggested to he U3ed to avoid the determinant singularity. 

The third one, are wings with both supersonic and subsonic leading 
edges, diaphragms have to be used for this kind of wings. 

If diaphragms are used, both values of the velocity poten- 
tial and its normal derivative on the diaphragm elements 
are unknown. However, two independent integral equations are 
obtained for each element of the diaphragm, one relates , and 
(df /?U ( ; to the upper geometry of the aircraft and whole diaphragm, 
the other one relates the same quantities to the lower geometry. 
Therefore, the total number of equations is equal to the total 
number of unknov/ns, then a system of algebraic equations can be 
solved for ; of wing element and and of the diaphragm 

■ k. i L> 0 

elements. If the problem is symmetric in the z-directios , then 
Pf/W is zero while k | t is unknown for each diaphragm clement. 

If the problem is antisymmetric in the z-directi^n, then „ . is 
zero while * , !)' is unknown. For these two cases only one inde- 

M> 

pendent integral equation can be written for each diaphragm 
clement, and the number of unknown of the diaphragm elements 



reduces to a half. 


1.7 Pressure and Generalized Forces 


The pressure on the surface of the body ia evaluated from 
the linearized Bernoulli theorem 




(1.32) 


Then is given by 


Ck . .P-P, . - .2 (?<g . ^ ?i ) 

r J e y.' J fc S’*' 

•S V** •» M 

For steady flow by using Eq. (1.4) Eq (1.33) reduces to 


1 P 


2 


(1.33) 


(1.34) 


(1.35) 


For oscillatory flow, ■ is derived in Section 4 as 

r 

4 < ’ ll . - * f0 + ;nk 4 > ) 

- X 

The evaluation of the generalized forces is considered as 
follows. The generalized force is defined as 

Q r - j'p i * 0^2 (1.36) 

where is the force acting on the surface of the body and 

*® ’ 

is the vibration mode. For lift, ' r l< and thus 

[• (J : - p n • k - . p 

Therefore 
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r -< .? p p (Ax 
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For pitch moment 
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since , \ . *a i r\ is negligible, therefore 

M * - jj' 0 *- x* ) d* d*\ 


( 1 . 40 ) 
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SECTION 2 

NUMERICAL FORMULATIONS 

2_. J Introd uction 

By introducing hypcrboloidal element, given by Eq. (1.27) , 
the integrals in Eqs. (1.30) and (1.31) are evaluated analyti- 
cally (See Ref. 3). In Subsections 2.2 and 2.3, solutions for 
b, and C are shown to be valid for any planar quadrilateral 
element inside the Mach forecone. For elements outside the 
Mach foreconc, b « C *■ ), since all these elements have no 

influence on tba element 7 . v . However, for elements inter- 

\ 

sected with Mach forecone, singularity problem arise , there- 
fore, s. lutions have to bo considered separately. In Section- 
solutions for this kind of element are considered. 


2.2 Source Inteoral 

. ■ ,i ■ »»— ■ ■ . - »■■■■ . . 

For element, inside the Mach forecone, the solution of b , 
in Eq. (1.31), is given as 


■- -i LV' 1 )- ?»<• ■-» » - 


( 2 . 1 ) 


whore 
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with *cc L^*(h'27)) 
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Note that Eq. (2.2) may be rewritten as 

i s « [- f <*.®" + 

_ / .7, ( W ' ) 

f* III <* !! • J. xAj./ (2.4) 




(2.4) 


where 


h * *>*« / |S, **t | 


* ' ' (2.5) 

is the unit normal. 

The following is to prove that Eq. (2.4) is valid for any 
quadrilateral planar element. For a planar element, the unit 
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normal n is 

independe*' 
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Furthermore, as shown in Appendix I) , (see Eqs. P.5, B.6 and B.7) 
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Uninq Eqs. (2.8), (2,10) and (2.12) one obtains 
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(2.15a) 


li<u iy, 

f *,M} = §■ Cf ^ 

• * d 

X '" ] h ||f if >ff/ If // 12. 15b) 


'/J ||f if If 

In addition as shown in Appendix C, 
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Finally, combining Eqs. (2.4), (2.14b), (2.15b) and (2,16), and 
noting that * yields 
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According to the Second Super-rule (Eq. A. 11 with a q, 

b « a, , c *= Q 2 ) 
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i' in 9 



I <7. > Q t \ 

lj« 

Substituting Eq. (2.19) into Eq. (1.52) and noting that H - 
for the element internal to the Mach forecone, yields 

i i ' i 

"'/lb r 
hfc 



( 2 . 20 ) 

1 


( 2 . 21 ) 


therefore, 

- n K* - M 1 *')- V 1 *-')- i|Ki) i,C-i, j) 

with I s given by Eq. (2.2). 


2.3 Doublet Integral 

Consider the doublet integral in Eq. (1,30) which can bo 
rewritten as 



(2.23) 


Note that !! ** 1 , for the clement internal to Mac i forccone and 
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Therefore 




( 2 , 24 ) 


tj7V c<> '' V 


Thus 


r — 1 A/* * 

((.q-) ,/( r 

* )j j X D * ' ^ " i D (k '" 1 ]* ■ L 0 C' *, 1 ) a 0 


(?. 25 ) 


( 2 . 26 ) 


whore 

9 *T 0 r . a ' * "*■ 
a5^1 Hffl 1 


( 2 . 27 ) 


In Appendix C, it is proved that 
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FINITE PART OF INTEGRALS 

3.1 Introduction 

« m ' r« • « r"*« ■' • ■* * t 

In order to extend the solution given in the previous sec- 
tion, the finite part of integrals are investigated in inis 
section ai.J the solutions for b and for the elements inter- 

sected with the Mach forocono are considered. In Subsection 

3.2, a general integral furction is considered. In Subsection 

3.3, the source integral is considered. In Subsection 3.4, the 
douhlet integral is considered. 


3 . 2 _ ) Gone r. ._1 Integral junction 

Consider the following integral equation 



(3.1) 


• <* 

where J (x) is a regular function and H(x) is the Heaviside 
functi on. 

It is observed that as x approaches zero, the integrand 
in Eq. (3.1) become singular. However it car. be proved as 
follow;} that this singularity is avoided after the integration 
is carried out. Let £ be an infinitesimal quantity, then 
Eq. (3.1) can be rewritten as 
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30 0 30) ^ ZO) 


vr v€ 


- TT 


C3.2) 


Therefore, the singularity contribution disappears and 
shou, not be taken into account. 


3.3 Finite Part of Source Integral 

By combining Eq. (2.17) with Eq. (2.20), one obtains 

Ifiiiil. -L (i . 1 - V 5 °' n . j „ *, o-k L?4f, ) 

ii tii n oK I mu a p 1 (i %) ‘ 

ai.O’K j °.±L ) 

<f®f /" ' 


( X q t o ft ^ y x P, on _ » K 

' H« H tii 


•+ t * ( v i^ii> 


(3.3) 

Therefore, the source integral, b , can be separated into three 

hk 

integrals, i.e. 

• i 


1 hk 


t r 


(-1 


-< -i 


h f)H 


I I 


. ~ + (3.4) 


where 


— JM± 


fl D xP, O'v - - - t 0a i 

1 - 0 r ^K^on h 


H i« 


5 2 = 
and -i 


a, x fti 

. « in 


O ll 


II $11 

LiA£l.j,a.oK 

(l$|l 


ll$ll‘ 


^(3.5) 


o a, 


Hii-j 


‘*T ^ (3. 


6) 
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S> r 


J'f «[« 

-I -I 


• TV 2 d ^ 


H 1 * J ' (3.7) 

Consider the first integral in Eq. (3.5) and assume that I,, 
is tho solution of S, through the double integrations. Then 

5, c T„ (1.0- I», Ci,h)- X„ C-«. 0 H,-«) (3.8) 

By using Eq. (2.15), Eq. (3.5) can be rewritten us 

*•- "i'f (3 - 9) 

Note that 'sj^0\ can be expressed as where 

yx 5. O h /^(>| ) is a regular function and is defined 

such that \ °% z O . Compared with Eq. (3.2), it can be 
concluded that after the first integration, solution along the 
intersection line of the element with the Mach Corecone yields 
no contribution, i.c. 

l„d.i) - o 


(3.10) 


or 


1„ (l.-O- I„ (-l.-l).O (3.11) 

if the edges of ^ s I , or ^ = - I is completely outside the 
Mach forecone. Otherwise, I s« is given by the first term of 
Eq. (2.2) if the corner point is inside the Mach forecone, or 


I„ U'V)- o 


a, o a, .■> o 


^ a, oa.<0. 


12) 


if tho corner point is outside the Mach forccone. 

Similarly, by using Eq. (2.14a), Eq. (3.6) can be rewritten 
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as 

51 ■ JA Jl.3? [ h * 5 * 051 } VPT 1 * f «•»> 

For the same reason given for S, . if is expressed as 

5 Z - Isa. C I » 0 “ Isi(»*-0- I** (-M)^ 1*2 f-^'O 0.14) 

the solution of S| is given as follows! 

In(l.O- Is*(l,-l) = O (3.15) 

or 

i« (-i.O- r« K-0- o «•«> 

if the edge of "f *= 1 * or f -1 .u completely outside the Mach 
foreccne. Otherwise, I*x is given by the second term of Eq. 
(2.2) if the corner point is inside the Mach forecone, or 



(T, <T l)= o 


■> a* Qd, ^ O 
i < o 


(3.17) 


if the corner point is outside the Mach forecone. 

e * 4 * 

Note that the values of ( and are evaluated such that 

^ 0 ^ - O . The solution of ^3 is considered together 

with the solution of c in the following subsection. 

hk 


3.4 Finite Part of the Doublet Integral 

The doublet integral is given by Eq. (2.26) as 

Uk = l[r B 6.i)-I l) (i.-i)-T,(-i.O-* T.f-I.-O] 

one can also express 5j as 


^3 = Is* 0 , l) - I51 ( I » “•) ~ Is ? (• I ) l) + It j (" I* - | ) 


(3. ID) 
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In Appendix C f it is proved that (Eq. C.15) 


ai„. 


- %0l l'*4.K O 


(3.20) 


lrt»*.ini{il 

Therefore, the doublet integral and Sj can always be expresscl as 

„. 21 , 

-l 1 — — r.l a A 


and 

5 


npiH 


• j \J ~ 


22 ) 


Roth integrands of the above equations are of the form of the 
integral in Eg. (3.1). For the same reason given above, one 
obtains 


Xll( It*) - XsiC‘ , iO r 

0 

(3.23) 


o 

(3.24) 

or 



In To (-<>-') 

^ 0 

(3.25) 

I,0*-l)-T, 

r O 

(3.26) 

if tho edge of ^ ? 1 or 

is completely outside 

of the 

Mach forocone. Otherwise 

, Ip and Ijj are given by Eq. 

(2.27) 

and tho third term of Eq. 

(2.2) if the corner point is inside 

the Mach forcconc, or 



^ i •* f($0«. Xjj OaJ({.fi)] 

(3.27) 

Id (f> V) - ** 2 . [ 

M ©*«X£ ©a-*) ( *■)] 

(3.28) 


As mentioned before, the value of is evaluated such that 

% 0% - O . 
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SFCTION 4 

SUPERSONIC OSCILLATORY FLOW 


4 . 1 Integra l Equation 

In this Section it is shown how the results obtained 
in the proceeding Sections can be extended to super- 

sonic oscillatory flow. Introducing the variables 


Xri V. i. Z: ?, T- Sl.»l 

B't i (4.1) 

and the complex potential £ such that 

4(XX 2)e 

(4.2) 

the integral equation for the subsonic oscillatory flow is 
given by 

4 = & fil- I MfaMiJt + 2.(± «{B M t))]wz 
™ JJ I art' ;'',i 1 sw e w;« ' /] 


(4.3) 


where Z surrounds body and wake. 


4 . 2 Round ary Condition 

The boundary condition is given by 


V 54 ’ 

V' 7' 




> t 

t. 



(4.4) 
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or 


, d rj A , * 4 J- 22 + AL l ?* 2+- * 0 

^ ^ r) * r £ >x J^ 1 <>* 




where Cj> and arc such that 


(4 .5) 


(f ; ^ * 0*** 


(4 .6) 


Ncx 1 ^ assume that the motion of the surface consists of small 
harmonic oscillations around a rest configuration, that is 

,vi r 


AK ^)- 1 Sl*'?' z ) e 


Then, setting 


4 >(W)* 

one obtains 


/JlT 


(4 .7) 


(4.8) 


.AT 



- ! ^ 

\ 

f a 

y \ 

,7 JIT 

- V Sol' , 

<? 

e 

'**>2 *rz 

* / 


+ lf?S, 4 

25 

,JLt 

•C 

B l 2X 




, ^ .723 e* 

/A 


fir 


« h* ' .■ 

T?/*x ax ax ax ax 2 a j 


o± 

ax x 




i2 a r 


1 c 


= o 


(4.9) 


I 

Assume that the surface is given in the form 



(Upper surface) 

-[z-z e ft*'|-z, (xyj^lo 

(Lower surface) 


(4.10) 

with 


!l • z H t i'x y j . out 

(4.11) 

1 ( (k?) . of*»; 

(4.12) 

or, in general by Eq. (E.7) with 


24 . 2Z„<. , c? & ; 

0X <^X 

(4.13) 

(4.14) 

-< 

>- 

it 

o 

(4.15) 

and 


2. -* (yy). oU", 

(4.16) 

M-- °o i ) 

2X 

(4.17) 


D 

D 

Q 

0 

I 

I 



Assume a ] bo 


SI . O(i) 


(4.18) 


This implies (see Eqs. E.25 and E.26) that 

<k ■ (4.19) 

Y ’ ° {d ) (4.20) 

iiJlT 

Neglecting the terms which contain t (which arc of order 

4 

f ) and separating the steady from the oscillatory terms, one 
obtains 

UyU*f ° f*- ZK 2X 

_ V? & & 'V Sg V cA t ~ *A$ + 2l 

+ /■(’ /.?,> 3'A , -oS of. ) _ o 
? ( ?X ax 3/ **/ ' 

Introducing / such th.it 


(4.21) 


( 4 . 22 ) 


S'* 


'} - • 7? (A / 

<p * </> e 


(4.23) 


Equation (4.22) reduces to 

/\ fl t\/ j S » • Jt 7-1/ r^t 

-V„AoV„ t 4> e . .-art /• « - 

I* 


.* iaS ‘ s % • & © (fl • 

] 


, f>J>. It- 7 . 0 

.:>/ ;.>x- ' 


(4.24) 
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Finally, neglecting terms of order i in Eq. (E.21) and terms 

3 


of order i in Eq. (E.24), one obtains 

-P 6 G V i r - 

W * <'t £i 

y . £ a, V <{> » - fl &S2- S * - JC 
- '***• Y *rt r \ jST IS J 


n 


(4.25) 


(4.26) 


In particular, for 


4 


' /** *•< Ni-KtW" 1 ] 


A .27) 


(where the upper ( lower ) sign holds on the upper [lower) surface), 
one obtains 


and 


✓A 


dN c 


£ * 1 j? [*• 

(4 .28) 

2 -- r j 

(« .29) 

1 r / A/,/ r i A/, 

IV g 1 2 I 1 

(4 .30) 

1 K/, j 


. E„4©p,,, <P __ a/ (■> k? + if. 

•MMX 
1 c 

• 1 r “> 

(4-31) 


whore 


{.'* *' Q 


/ 


A 1 


- id 


( 




(/. 3?) 
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4.3 Pressure Coefficient 

* ■ ■ 1 ■ ■ < « — ■■ M I 

The pressure coefficient can bo evaluated by using the 


linearized Pernoulli theorem, as 

v-|. 

- 

v £/ ) 

" *>x / 

. -2 

f l 

, / 


» 7 * 

£5 ^ X 

For occillatory 

flow, setting 


a/ : 7 r ifil/T-rf*) 

<jj * <p C ' 


A-> 


;*/ir 


T 


(4.33) 


(4.34) 


(4.35) 


one obtains 
rJ 


AJ 


<r 


3/ \ 

^ ) 


- 2 (i ■»* J i 

.2 fo/l - * ) i V 7« 

L ( /-I ; d> ^X J 

- /x 1 . 72 ** 

- -2L 1 * t T 5 '- / € 

£ / a/ ✓'X J 


. J2 



(4. 3C) 
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SECTION 5 
NUMERICAL RESULTS 

5.1 Introduction 

The formulation presented in the previous sections was 
imbedded into a computer program, called SOSSA ACTS (Steady 
and Oscillatory, Subsonic and Supersonic Aerodynamic for 
Aerospace Complex Transportation Systems). Typical results 
of this computer program are presented in this section. 

5.2 Rectangular Wing in Both Steady and Oscillatory 1'low 
The results in Figs. 2 to 4 are relative to a rectangular 

wing with aspect ratio AK «* 3 and with a biconvex circular arc 
section, 51 thickness, with sharp loading and trailing edge;. 

Fig. 2 shows the distribution of the pressure coefficient C p 
on the lower and upper surfaces of the wing with ^ ■ 0 4 and 
M «= 1.3. Fig 3a shows the distribution of the lift coefficient 
on the wing with cs = 5° and M ■ 1.3. while Fig. 3b shows the 
distribution of C^ on the lower and upper surfaces of the wing 
with e\ «=> 5° and M = 1.3; these results are obtained with 
NX =» NY « 7. Fig. 4 shows the distributions of the absolute 
values and phase angles of the lift coefficient C . of the same 
wing oscillating in bending mode 

l * 0 , <6043 - / 7 o?S£ I >•///. I* - • 1**6* I't/bl'-* ■ ?*** 7 I 

with K « 0 * /$ u» - .1, M - 1.3 and NX « NY = 10. All the 
above results are compared with the ones obtained by Lessing, 
Troutman and Menees (Ref. 4). The results obtained for Fig. 2 
are also compared with the analytical two dimensional solution 
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which can bo easily evaluated. For, this problem can be 
treated as a two dimensional problem in the central region 
of the wing. For this case the pressure coefficient is given 
by (Ref. 4) 

C y r -ll?- (5.1) 

For biconvex circular-arc wing, the equation of wing section 
is approximately given as 

? s- -?t (X - P >*) , ( - o.f £ X £ O. ) (5.2) 


Noting that T = 0.05 and 3 * 0.83, one obtains 



7 


- & X 


i.e. c varies from -0.24 to 0.24 linearly. The results 
P 

shown in Fig. 2 are in excellent agreement with Eq. 5.3. 


(5.3) 


5,3 Convergence Analysis 

The convergence analysis of the problem considered in 
Figs. 2 and 4, is presented here. The distribution of the 
velocity potential along % r o for the problem for Fig. 3 
(for different numbers of elements) is shown in Fig. 5. The 
curves are obtained with NX = NY = 5,6 and 7. These curves, 

indicate that the convergence is very fast and that 144 
elements on the whole wing, or NX ■ NY = C are sufficient 
for an accurate analysis. For oscillatory flows, the distri- 
butions of the real and imaginary parts of the velocity po- 
tential along 2 y/j, to.JT for the problem for Fig. 4 are shown 
in Fig. 6. The curves are obtained with NX = NY 6* 
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nnd 7. From these curves, it is safe to say that 144 elements 
(i.r NX - NY ■ 6) are sufficient for an accurate analysis. 


5.4 Delta Wing with Supersonic Leading Fdge 

Fig. 7 shows the distribution of lift coefficient per 
unit angle of attack for a delta wing with supersonic leading 
edge and 

lil s & / A r I. 2 


where A is the sweep angle of the leading edge. The results 

obtained with NX = 8, NY «* 12 and M ■ 1.2 are compared with 

the exact solution which is given by (Ref. 5). 

f A ik »•' n r ' l- *i f i+mt 

** C P ■ , T ; y-~s=s R> Coz ----- - tre 


ih « 


~ ] 
r J 


where f ■ £>[ . The numerical results obtained are remarkably 


accurate. 


5.5 Wing-body Confinu. . at io n 

The present method m general enough to extend to any 
arbitrary configuration. Following is an example of this 
application. A wing-body combination in supersonic flow with 
M = 1.48 is considered in Figs. (8a) and (8b). The combin- 
ation is composed of a wing with chord C « 3, span S= 9, 
thickness 7 = 0.05, a forebody of length LA = 6.0 and radius 
varying from 0.0 to 0.75 linearly and a midsection of length 
L m = 3.0 and radius r = 0.75. Wake and at 'body are not con- 
sidered. The angle of attack of the wing is = 1.92 , 
while the angle of attack of the body is o\ p = 0 , To obtain 
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the results, 580 elements on the whole configuration (NX - NY 
10 on the wing, NX * 5, NY « 3 on the body, NX * 10, NY ■ 3 
on the middle section) are used. In Fig. (8a), the distribu- 
tions of the lift coefficient per unit wing angle of attack 
along chordwise direction are presented. The curves are 
plotted at different values cf y/r, and are compared with the 
experimental, as well as analytical, ones obtained by Nielsen 
(Ref. 6) and Woodware, Tincoco and Larsen (Ref. 7). In Fig. 
8b, the distributions of the same quantity along fuselage 
at different meridian angles are shown. 


5.6 Com puter Time 

All the above results are obtained on the IBM 370/145 

available at the Boston University Computer Center. The 

computer time for the problem for Fjgs. 5 and 6 are given in 

Tables 1 and 2 respectively, where N,p is the total number of 

the elements on the whole wing. N t ^_ is the number of e- 

quations to be solved (by using IBM Subroutine GELG) , and N, 

is the number of diaphragm element. 

TABLE 1 Steady Case 


n t 

NX, NY 

Neg . 

N d 

Computing Time 

64 

4 

28 

12 

11 sec. 

100 

5 

40 

15 

22 . 8sec. 

144 

6 

54 

18 

42 . 5ccc. 

196 

7 

70 

21 

140 . 5sec. 

— — — 
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TABLE 2 Oscillatory Case 


n t 

MX, NY 

Neg. 

N d 

Computing Time 

64 

4 

25 

9 

21.4 sec. 

100 

5 

34 

9 

38.0 sec. 

144 

6 

45 

9 

65.4 sec. 

196 

7 

58 

9 

118.9 sec. 


It should be mentioned that the advantage of symmetry with 
respect to z and y was taken. Therefore, the number of unknowns, 
or equations, to be solved is only one-fourth of the number of 
the elements on the whole wing plus the number of diaphragm 


elements . 
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Figure ^.Analysis of Convergence: Distribution of f = ?e'‘' ,v,/x Versus x/c, at 
?y/b “ 0.5, for Rectangular Wing With Biconvex Section, Oscillating 
in Bending Mode in Supersonic Flow, for AR = 3, t c 0.01, M = 1.3, 
K * 0.1. L *-* 0°. Nrv = 3N . 
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APPKHDIX A 
SUPFRALGLBRA 


A. 1 Super-product 

As mentioned in Section 1 , in order to simplify the algebraic 
manipulation for the supersonic flow theory, it is convenient to 
introduce a special algebra, called supersonic vector algebra 
or su per- algebra . In addition to the rules of the ordinary 
vector algebra, the super-algebra includes a supersonic dot 
product or supe r-pro duc t 

Cl © B » A* bjr « A) t>g . 


The additive and distributive rules are obviously valid for 
the super-product. Note that a o a is 


ft O 5 ^ o Ur 0 K J (A * 2) 

tha< is for a pointed, respectively, inside, on, outside the 
Mach cone (Fig. 1). Hence, in addition to the ordinary norm 
of it vector (or dot-norm ) 

l£ I s J a • 0 (A. 3) 

it is convenient to introduce the supersonic norm (or super- 
nqrrn) 

a (1 * 



Fin«ally, it is convenient to introduce the concept of covcctor 



(A.f») 
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With these notations, it. is immediately verified that 

a o 5 * a c . h s a • b c (a. 6) 

It may be v/orth noting that 

Q . IT * C « C fc O h * C * a o b * cf S ft o 0 C * C < (A. 6 a) 


A. 2 First Euper-r u le 

Throughout the subsonic finite-element formulation (Ref. 6) 
the following rule is used 

^ * p)» fc.» (' ) * (&* £)(!?• 3) - (3 • d) (b ’ c) (A. 7) 

The corresponding supersonic rule, called for convenience, 
first super-rule , is also valid 

(a>t)o(e,J)' (ao c)(Eeo).fao 2 )(ioc) (A 8) 

For 

( d < t) o (c * 3) 

=• ( Gy bg - 0 S b y) ( Cy d s • C t dy) 

~ ( Gt bn- Qx b* ./ ( Ce dx - Cx dg) 

~~ (Qx by - b*) ( C* Ay -Cy dx) 

~ bt Cydf Cl J d/ ** !>/ Cf di 

by Cg d a — G * bg C r d g Clg \>h Cm di -t Ox bg Ci dx 

~ 0 K b; C A dy - 0, bx C 1 d* r h C 1 ^ ^ (A * 9) 
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whilc 


( a o t)( l> o 'd)-(ao djfbo c ) 


— (a, Cm - flf Cf - fle C i)( l j, d A - by dy - b{ d*} 

- Ay - Q* ^i)( b* c r - b /0 - £*) 


r= Q Cm ( l>/d* - bydy - I'rdz) 


/ 


- f bx d* b*dz) 

- A*Cz ( b/CL - 

- d, (W* - JyO- ^ 


+ ( 3 , d, (t, C, -M>- brC') 

* fig di ( bs Cm ~ by Cy - biC^) 


tj Cy di + fit by Cidf- fly beC* Ay - by C f &i 

— i>x d* - hi Ce d s *■ fl e bjt Cf d^ * fl* be ft d* 

- by Cy Ay - fly b t Cy A x ■* Qgby CjAja* Ay b* Ay Ay 


(A. 10) 
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A. 3 Secon d Super-rule 

A second rule of the super-algebra is 

(ao 0t)( l » c) O (S< c) - ( a ■ l* o ye) 

sr a o r (It ~c) ° <t) + ao L(c<Z)o(cxd) 

^ otoc (G*c)o (StS) + J Q I 

(A. 11) 

Note that the dot product appears in the triple product. In 
order to prove Eq. (A. 11), consider the regular vector algebra 
rule 

0 t (l * C-) s: l ( 4 » c) — C ( G • T>) 


This yields, for the covector a e , (see Eq. A. 6) 

< (l ) * C ) *=■ F f tf c * c)— C. ( G • F ) 

* F (go c)~c (a ot>) 


(A. 12) 


On the other hand, according to Eqs. (A. 6) and (A. 8) 

Q c , ( If l) O ~d* K ( l f C ) 

B- a c o C£*c)o ( d‘o (l* S>&sc)j 

st ci o d ( F* c.) o j . , 

while, according to Eqs. (A. 6) and (A. 8) 

( i ( d o c)- c (5 o (I (doc)- c Cato'S)) 

- b(-)L> ,\<z> c a&C — A G> C d<£> >» 

" * © S> tjji C -V- 3«C i\o o ,V,>© 

~ 3oc (1 >mCoIm a) + aoT>(c,loc, a) 

t O C !•> a c r »X 4 * .( • u C. a £> 0 Q * ■ I 


(A. 14) 


1 
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1 

Combining "qs. (A. 12), (A. 13) and (A. 1*1) yields 


I 

21 o A b x c J 


1 

T 

= ( (t*a) o C Jc, ( i , £)) 


1 

- (Hoc) - c 


•1 

b flee f ‘5 » ce\B * dl) 6o5 ( c x b]oS k a!) 


D . 

^ * 

z 4 a> c ( o k cjo! 'o * a. j + 5 © I ( c v i>)o(c‘/r.i j 

(A. 15) 

ii 

that is the second super-rule, bq. (A. 11). in particular 

$ 

u 

for aji/iq, T> s a^, c 5 a^, one obtains 


0 

p j- a,* \ - 1 f- a,* Q t 1 * 



* 1 S,'0* o *,'j.+ l*£, (a t 'Q, 0 l' j) 

(A. 16) 


A. 4 Third Supersonic Rule 



A third useful formula, called the third supersonic 

rule , 


a * <} 0 b * c a • { < y + a* ■£ 0 b * c a • j * 6 


Li 

(k * 0 b f C A * 6 • -f - 0 

(A. 17) 

w 

The proof of this rule follows 


D 

( a * 3 0 l < c)(a • { * <j )-( a > j 0 1 * c ) (a • f • 3 ) 


II 

- ( a oJ> c 0 3 - a 0 c lo 2 J a * ? • j 


IJ 

n 

- (flot coj - a OC lot) r f • J 


11 

TI 

r a 0 1 ( c 0 2 coy atj. 2 ) 


11 

~ a 0 c ( l 0 i a * -f -j - 1 0 $ a » 3) 



1 

I 
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- a • V (c*. 6 a » (.y .. c *• f a* f -d) 

- ci ' C' ( b e > d * -f • $ — b c - <j a > { . J ) 

= a-V[c c * ( a * f ))• 9 * ~J - a c‘ []> c t ( a 0 f )}• 2 'j 

* a- ]> l [ (c«- d» (?'•*)( f • d* f ) J 
e [ a • j> e c *• -f - £ • c * ~y L > ) a • 3 * 

e (aoj/cof^a o c b o f ) « ( Q • { o b * c ) Q •(!*<) 

(A. 18) 


In particular, for a ■ b - q , c « f «= a 2 , d «=■ p 3 , g *= a , 
the third suporrule reduces to 


(jM® s)- 

-d‘ i) 


(A. 19) 


A. 5 Fourth S upersonic Rule 

In subsonic theory, it is easy to show that 

I"' j *4/ * * // xi/ If 


(A. 20) 


The corresponding equation for supersonic theory is called fourth 
supersonic rule 

+ ( f. rf i 0 p4 ^ I If ' # H 


(A. 21) 


The proof is shown as following, according to the second 
super -rule. 

i© i (f- ) * (f< *!<?}*«•) 

r | o % C ( ° i & *~ Q ‘ 0 ^ ^ ~ 1 ° *' * *' ° *' ' ¥ 

a 

- f 04, 4, * 4, 0~°> *1 0 + ('1° t *• ° ** " f° 4 ^ 

= (<j.0$) [c{sQ Q, 4f. o\ f ] -'J-O f)' ? l 0 * - 

[ , Q t G\~4,0 j! 4, (>'<?,]- £ 6» j j o 4, (Ci04x 4,C-£- 

+ r f ,p ®*>'?** *?«*•> 

t IX q\ J # 4* ) 

- [ | ^ £ Q,0G, - { f 6 4t ) ][ fof 4, 0 4, ~(j 04,) ] 


' J Oft ‘i’i ) 
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APPENDIX B 

T> f i > F> 

— A HP — — 

V 


In order to prove Eq. (2.1) is valid for quadrilateral planar 
element, it is necessary to p*-o/c the following two equations 


■> f, i 
* f = ii i a 
if. _ i 

*7 " iif ii 

Proof of Eq. (B.2) is shown below. Proof of Eq. 
except for the fact that a^, is replaced by a^. 
different cases in t' 2 . Consider the first case 
Note that 


(B.l) 

(B.2) 

(B.l) is similar 
There arc three 

G > O 


since 

and 


14, H = o 

c b‘fh 
if } 

, , (d ,z)tKfj,«) 

- * ( 4 . < 4,)0 (j 


o 





(B.3) 


(B.4) 



= fltir f 

I I /ioi /— r - -\ 

* »<,i Ifli&ii+jo* VT^r^ C, ‘ 4 ^ ^ 

a ?,(! a fnt,» ~}c~Z C~(l° 

I 

Consider the second case, Q t O Q, = o 

T> F» _ _J>_ ( ftjj 1 

Z >-9 “ D/; l ^ G> 5 , J 

4 J i 3 U 

= — 1 1 II 

^ O Qx *') 1 

i___ t o i _ J 

j °*> , rpT " (, f* 


since 


53 


Consider the third caso, Clj C < Q 



J I f » <?. H {_ - a 

11 J * 11 ' 2 " 

ti 4, 1» 

J('f o \i(r4*o*») 

I 

»?• 

Combining (B.5), (B.6), and (B.7), yields 


»_F> _ t 

<"] ~ ijii 

Similarly, yields 

& F> _ J_ 

* Ilf If 


<?, o £ 0 


Qy O 4y 


£ 


(B.7) 


(B.8) 


0 


(B. 9) 
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APPENDIX C 

DERIVATIVES OF Ip AND I g3 


C . 1 Introduction 

In this Appendix, it will be shown that the second mixed 
derivatives of 


If • nl If I l| ■/),. a, I 


1(> '> ftfll 


are given by 


and 


T> J, 


\ • ax* a> 
■ tjt* 


S3 



•/? 

the third term 

v f? i 

-)f?, 



9 a. 

b 

M 

•t 




S 


Note that 


s l? 


r< 


(C.l) 


(C . 2 ) 


(C . 3) 


(C . 4 ) 


(C . 5 ) 


(C . 6 ) 
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and 


jlL * 

< 0 ^ 1 

(C.7) 

4 . 5 , 

(C.fi) 


C.2 Derivative of I, 


Consider Eq. (C.l) or 


Ip - ^ 


- 1 1 -%t 4, • 

if ii $; 


where 


- - , = '>>/* ( j-n ) 

If • n 1 ' r 

The derivative of I D with respect to is given by 

lli -2- c tfl (s3.< *'.?]*«' 

*1 *1 


* - 5 ; 


ffc ] C ( fl < l ) 

4 

7771 - o/ ; \ y 

' Ji e i 


(C . 9) 


(C. 10) 


I [ #i< 4, ° P + i < fro f / 4, + J.y 4iO a, >: 4 , ] ~r — r — - — — 

/ 0 J 0 

f * q , 0 i* q> r\o~i - - - r- — r-v- ----- 


X 
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£ o j ( j • 4, k 4 J 4 

(|#| (p f if 9f9*(f*frlt 

X J ( a, < ^ o Q» + / * fi o 'f X ti>) 'f O Jf / • Si * 4 

- £<«•*/ *4 [ «>oj ;.*<«, +joj ;. p*Z t ] 




(?0 j Uf *4, * Q, )'+(}* lop * yjjjf 


•+ [(?' ft oj<i)pi t <i ~(j**iO j,i)f 



(C.ll) 


Next note, ”8 shown in Appendix A, Eq. (A. 20), 


1'J <c. 


12 ) 


Moreover, note that (see Eq. A. 19) 


( ffl 4 vtf, o 

ph y{°i’ ( p 

<?, O P 4* 

) 0,0 Jfjj ■?, 

-r <? 

o r j roM) **?/ )~ 

■ (p«< c>~( 

- «j)(fp, ‘ 

VJ t 

= f7*>ot 

a, 'V, - 4* O Q> 4 


f 


-( jejf / C. - ] o'a, | 

0 3. ) 4, < 

>? /?•* 

K 


fi 0 Vi, - £ 0 f 

op,) j- 

4“i *■ <?i 

~) 


-(in 

q x <>4, / 

C 4i ) -j - 

b **jf 


• ~j t°i 

4, o - C? <> 5, / 

) fef, 

f • *. X 4i. 



(Q\ (. ^ f 1 } ^ 

* At O 4i j 

' b * 4* ) 

-) 


jQOi 

(fop, }■«<<*. 

no o-i 

)j 

fop 

m -\j°i « 

'.oo, - if oa.fj 

I** }■ 

a 4x 


+ (]>«. 

o f- o.) if-#,* fil'f o ] ) 



= if n, n 

J ( jo 4, % • 4, A \ 

-M i 

• 4* x /' ) 



(C . 13) 


since 


0 f Q t 04>-('lG4,)-.j*4 i G>^4 i - 'Up <?, if 


(C . 14 ) 


Finally, combining Eqs . (C.ll), (C.12), and (C.13) yields 

dI „ Ilf pr u <r^ .J'ff.i ird 

7 >') Ilf- 4, \\‘ Ilf • C. n‘ II; II 
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^* 4,1 IjJ 


?, joj j' 4 "< ) 


(C.15) 


Next, consider the second mixed derivative, noting that 


-p ('/<§,) - s -y[ if, + ')[>)] - o 


one obtains 


3 I p 

?f’7 


_i _ 2 _ f Jof I- to« i ■*<*«> 

tj'ir *f \ j Jfoj 


■s? —-.r-r, ) - 1 ** ^ C - i'C\y }“!<<<?>) 

11*4X1 (^V 1 t 


~r 


==r £o <yyS4*‘fi t ^ * ft - 4 e> 'u *} 


Fjjj • 'sr ' ' ■ ■ N 

- £ <? <?, ' ?» - % ®lKl ■&•?>) 
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